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Electronic states of an Cg electron are calculated in the system composed of two MnOe octahedra 
with the inclusion of the Berry phase acquired by parallel transport. Based on this calculation, a 
comment is made on the controversy between "Wigner-crystal" and "paired-stripe" models for the 
the insulating charge-ordered manganese oxides. 
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I. INTRODUCTION 



In quantum mechanics the Hamiltonian H alone does 
not determine energy eigenstates of the system. We need 
information on the "wavefunction ip "which includes the 
boundary condition in general and for the case of a many- 
electron system, antisymmetry due to the Fermi-Dirac 
statistics. The antisymmetry may be regarded as addi- 
tion of the phase factor e*'^ to (p in the interchange of 
any pair of electrons of the same spin. Note that this 
factor e" brings about an important concept of "the ex- 
change interaction" J "which is just the energy difference 
bet"ween singlet {S = 0) and triplet {S — 1) states in an 
interacting t"wo-electron system. 

Addition of the phase factor to (p also occurs in the 
course of Berry-phase connection; the Berry phase is a 
geometric phase in ip acquired by parallel transport. R] 
On the analogy of J, we propose a ne"w concept of "the 
geometric energy" by considering an energy difference be- 
t"ween the states characterized by different winding num- 
bers w (the Chern integers) associated "with the topolog- 
ical invariant in parallel transport. 

In order to substantiate the above proposal, we calcu- 
late the geometric energy in this paper in the simplest 
possible system, namely, a system consisting of a sin- 
gle electron moving back and forth bet"ween t"wo Jahn- 
Teller(JT) centers. In Sec. 2 we specify our system pre- 
cisely by writing H and a prescription for parallel trans- 
port. The lowest-energy state for each w is given in Sec. 
3 and Sec. 4 deals with our perspective into the stripe 
structures observed in Lai_2.Ca2,Mn03 with x > 0.5. 0] 



II. TWO-JAHN-TELLER-CENTER SYSTEM AND 
PARALLEL TRANSPORT 



Consider an Cg electron in the system composed of two 
Mn^+Og octahedra with one oxygen ion common to both 
octahedra to provide a transfer path with t the transfer 
integral for the Cg electron in z-direction. We assume 
that the spins of t2g electrons are aligned in one direc- 
tion. Because of the strong Hund's rule coupling, spin of 
the eg electron also aligns in the same direction. In this 
sense, there are no spin degrees of freedom and we may 
write H in second quantization as 



H = -t(b+b2 + b+h) + E,T Y. [2d+(gf r, + 9fV,)d, 

(1) 



i=l:2 



,(2)^ 



-,(3)^ 



with the Pauli matrices, t^ and r^, and d^ = {a'J,b'J), 
where aj and b'j are, respectively, creation operators 
for electrons in dx^-y2 and d^^'^-r^ orbitals at Mn site 



,(2) 



.(3) 



j{= 1 or 2), q- ' and q, are dimensionless vibrational 
variables representing, respectively, {5Xj — 5Yj)/\/2 and 
{26 Zj — 6Xj — (5Yj)/V6 modes around site j, and i?jT is 
the static JT energy. Since wc shall treat static distor- 
tions in the adiabatic approximation, the kinetic-energy 
term of ions is not considered in H . 

The Berry-phase effect at site j is best included 
by writing vibrational modes in polar coordinates as 



(2) 



qj sin 6j and q 



(3) 



qj cos 6j . Using 6j , we 



transform the operators a, and bj into hj and bj 
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as dj = e**j/2[ajCOs(6'j/2) + 6jSin(6'j/2)] and bj = 
Q^Sjf^^—aj sm{9j /2)+bj cos{6j/2)]. In this representation, 
the second term (the site-energy term) in (||) is diago- 
nahzed. The molecular Aharonov-Bohm effect manifests 
itself in the phase factor e*^^/^ at each site j at which 
topology of the singularity is understood as a double- 
plane feature. 

The transfer between the two sites connects 9i with 
02 . We make this connection by parallel transport, by 
which we mean that if we start with a state represent- 
ing an electron at site 1, ipi, we make the state evolve 
from ipi in the direction orthogonal to it, namely, along 
the direction of a state representing an electron at site 
2, ip2- This parallel transport indicates that an eg elec- 
tron picks up only a constant phase ^o by a trip from 
site 1 to site 2, namely, 02 = 0i + Oq. A further trip in 
the same direction from site 2 results in the total phase 
as 6*2 + 6*0 = 6*1 + 26*0. Now in general, a two-site system 
is equivalent to a corresponding periodic lattice system 
with periodicity two. Thus the phase 9i + 29o should 
represent the same physics at site 1 at which the phase 
is 01 . This reasoning leads us to the conclusion that 29o 
should be equal to 27rz/; with an integer w. If we make a 
more general argument on topology, it turns out that 
w is nothing but the winding number, a topologically 
conserved quantity which is useful to specify an eigen- 
state of the system. Note that w describes a geometric 
structure as to how the double-plane structure at each 
site is interconnected. 



III. GEOMETRIC ENERGY 



Let us calculate the lowest single-electron energy for 





each w{= or 1), Eq'. We adopt a variational proce- 
dure by considering a wavefunction |^q 



,(-)\ 



I* 



M\ 






e^^^/2(a("')a++/3f)6+)|vac) 



\qi = 91,92 = 92; 



(2) 



where a^™ and /3,- are real numbers, |vac) is the vac- 
uum for electron operators, and \qi — 91,92 — 92) is 
the phonon wavefunction representing the (5-function-like 
distribution of qj around qj. We have optimized all 
the parameters involved in the problem except for 62 
(which is related to 61 through 61 + ttw) in order to 
make (^q |iJ|^Q ) minimum under the condition of 
(^[,'"^1*^'"^) = 1 for given t and E^jt. We have found 
that the lowest energy is obtained at 6*1 = irrespective 
of w. This is the condition for the e^-electron orbital 
to polarize in the transfer direction. (If we start with a 
different configuration, say two sites in x-direction, we 
obtain completely the same physical results with the in- 
terchange of the role of z-axis with that of a;.) However, 



l^g ) is different for different w in an interesting way. 
For w — 0, we obtain 



,(0)v 



i*r>-(M"^&^+/3rfej)ivac) 



:>(0)h+^ 



■Mi-^r 



,(0) 



>92 



-^r 



(3) 



with the coefficients /3 shown in Fig. |l|(a) as a function 
of EjT- As long as Ejt is smaller than t, symmetry exists 
between the sites, but for larger E^t the site-symmetry 
is broken; an Cg electron tends to localize in site 1. (Of 
course, the state with the interchange of the sites is pos- 
sible, but there is no overlap between these two states 
because of the orthogonality of phonon wavefunctions.) 
For w = 1. on the other hand, we obtain 



I* 



(1)\ - fRWh+ 



W, 






(4) 



with the coefficients /?} and aj shown in Fig. |1| 
In this case, the site-symmetry is broken even for an in- 
finitesimally small positive value for _Ejt due to the fact 
that the lattice at site 2 never deforms, because with this 
winding number, an eg electron at site 1 in the lower 
adiabatic potential plane is parallel-transported to the 
upper adiabatic potential plane at site 2. In this sense, 
w — 1 describes "the inter-potential-plane connection" , 
while w ~ "the intra-potential-plane connection" . 

We plot the corresponding energies -Eg as a function 
of EjT in Fig. in which we give an analytic expression 



p(0) 



p(i) 



for Eq . (We can give Eq only numerically.) The term 
— i^/2i?jT clearly indicates the hopping nature of an elec- 
tron localized by the JT stabilization energy — Ejx for 
EjT > t. We find that Eq is always lower than Eq '. 



rW 



(0) 



The geometric energy, Eq — Eq , is negligibly small 
in both large- and small- Ejt regions, but it becomes as 
large as O.lt for Ejt ~ t. 

Quite a similar situation occurs in the J-problem in a 
two-electron system in which the conserved quantity is 
the total spin S; the energy of the singlet state Eq^'^ is 
always lower than that of the triplet one Eq^^ . Note, 
however, that application of external magnetic fields can 
compensate the energy difference to stabilize the triplet 
state. Analogously, we may stabilize l^g ) by compen- 
sating the geometric energy with application of external 
stresses to the system in z-direction, because the state 
l^g ) is less distorted in the direction than j'l'g ). 



IV. DISCUSSION ON PAIRED-STRIPE 
STRUCTURES 

Considerable accumulation of experimental data has 
been made as for the crystal structures of the insulat- 
ing Lai_2:Caa;Mn03 with x > 0.5. Everyone agrees the 



existence of the charge and orbital ordering in these com- 
pounds, but a controversy continues about their detailed 
crystallographic superstructures; some claim the "paired- 
stripe" model IQ] and others the "Wigner-crystal" model. 

i 

In order to shed light on the argument, we have ex- 
tended our calculation from the two-JT-center problem 
to a two-dimensional lattice system, [pi The key idea is 
to identify a zigzag conducting path of Cg electrons along 
which we define the winding number w. We have found 
that the state corresponding to the Wigner-crystal model 
provides the lower energy than the paired-stripe state, 
though the energy difference is quite small. This result 
is obtained without any consideration of the long-range 
Coulomb interaction V. We have included the effect of V 
recently and found that further stabilization is obtained 
for the Wigner-crystal state over the paired-stripe one. 
Thus, in terms of energies, it seems to be certain that 
the former state is favorable. 

However, we have investigated the electronic state of eg 
electrons in the observed paired-stripe structure, evalu- 
ated the corresponding winding number w, and found an 
interesting relation between w and the Cg-electron den- 
sity specified by x as S 



1 



Number of Mn ions 
Number of Mn ^ ions 



(5) 



Since w is an integer, some specific values exist for 
x[= w/{l + w)] such as 1/2, 2/3, 3/4, etc. This reminds 
us of the experimentally observed special values of x and 
the related lever rule. M Thus, if the crystal is placed 



in an "environment" to favor the state with this wind- 
ing number w, we can expect that the paired-stripe state 
may be observed. At present, we do not know precisely 
how we can define the "environment" and this may be an 
important issue in the future not only from a viewpoint 
to clarify the crystal structures of Lai_a;Caa;Mn03 but 
also from a fundamental aspect as to the identification of 
a physical variable to control the winding number or the 
geometric energy. 
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FIG. 1. Coefficients to specify the lowest single-electron wavefunctions plotted as a function of Ejt for (a.) w = and (b) 
w = 1. 
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FIG. 2. Lowest single-electron energies as a function of E}t for w = and 1. 



